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Abstract 



This paper studies the perturbations of the continuously self-similar criti- 
cal solution of the gravitational collapse of a massless scalar field (Roberts 
solution). The perturbation equations are derived and solved exactly. The 
perturbation spectrum is found to be not discrete, but occupying continuous 
region of the complex plane. The renormalization group calculation gives the 
value of the mass-scaling exponent (3 = 1. 
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I. INTRODUCTION 



Numerical calculation of spherically symmetric gravitational collapse of a massless scalar 
field by Choptuik and subsequent results for different matter models and symmetries 
(see for example Refs. pj-^) spectacularly demonstrate that critical phenomena occur in 
the gravitational collapse, and that the near-critical behavior is universal in some important 
aspects. 

There are two possible late-time outcomes of gravitational collapse, distinguished by 
whether or not a black hole is formed in the course of evolution. Which end state is realized 
depends on initial conditions, specified by a control parameter p that characterizes the grav- 
itational interaction strength in the ensuing evolution. For p < p* (subcritical solution) the 
gravitational field is too weak to form a black hole, while for p > p* (supercritical solution) 
a black hole is produced. The critical solution, corresponding to a control parameter value 
p* just at the threshold of black hole formation, acts as an intermediate attractor for nearby 
initial conditions, and has an additional symmetry: discrete or continuous self-similarity, 
also referred to as echoing. For supercritical evolution there is usually no mass gap, and the 
mass of the resultant black hole is described by the mass-scaling relation 

Mbh(p) cx |p - pY, (1) 

where f3 is mass-scaling exponent. 

General gravitational collapse is a very difficult problem, even if restricted to the case of 
spherical symmetry. As a result, most studies of critical behavior resort to numerical calcu- 
lations at some stage, with few fully analytical solutions known. Finding critical solutions is 
made easier by the fact that they have additional symmetry; requiring solution to be contin- 
uously self-similar can simplify the problem enough to make it solvable for simpler matter 
models. Both the stability of the critical solution and the exact mass-scaling exponent can 
then be determined by linear perturbation analysis, as suggested by Evans and Coleman Q 
and carried out in Refs. [^|-|^ . 

In this paper we consider an analytical continuously self-similar solution of gravitational 
collapse of a massless scalar field constructed by Roberts 0. Originally it was intended as 
a counter-example to cosmic censorship conjecture, but it was later rediscovered in context 



of the critical gravitational collapse by Brady [10| and Oshiro, Nakamura and Tomimatsu 



rT| . Although this solution does exhibit critical behavior, there are two problems. First, as 
any self-similar solution, it is not asymptotically flat, and the mass of the black hole grows 
infinitely. This can be prevented if we cut the self-similar solution away from the black hole 
and glue it to some other solution to make the spacetime asymptotically flat and of finite 



mass, which was in fact done in Ref. |12[. A second, and bigger, problem is that growing 
non-self-similar modes will dominate the mass-scaling exponent if they exist, but they are 
not accounted for by the self-similar solution. 

This paper addresses the second problem by doing linear perturbation analysis of the 
Roberts solution. We calculate the mass-scaling exponent (3 and analyze the stability of 
the critical solution. The remarkable feature of this model is that it allows exact analytical 
treatment. 
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II. SELF-SIMILAR SOLUTION 



We begin by presenting a continuously self-similar spherically symmetric solution of the 
gravitational collapse of a massless scalar field (Roberts solution). It is most easily derived 
in null coordinates. Therefore we write the metric as 



2e^^ dudv + dVl' 



(2) 



where S and R are functions of both u and v. The coordinates u, v can be reparametrized 
as u{u), v{v), so coordinate choice freedom is given by two functions of one variable each; it 
can be fixed by setting E = at two null hypersurfaces v = const and u = const. With this 
choice of metric the Einstein-scalar field equations become 



(3a) 



9y^ 

R. 



R,uu '^'^,uR,u 



R{^,vy 

R{^,u)' 

R,v^,u , R,u^,v 



R 



+ 



R 



(3b) 

(3c) 
(3d) 

(3e) 



The first four equations (|3a| - |35D are the Einstein equation for metric (0), and equation 
is the scalar wave equation = 0. Equations ( |3^J3HD are constraints. Note that in the 
case of a massless scalar field the wave equation follows from the Einstein equations, so one 
of the equations ( Pa|J5b| , |5^ ) is redundant. 

We turn on the influx of the scalar field at the advanced time f = 0, so that the spacetime 
is Minkowskian to the past of this surface. The initial conditions for system (|^) are specified 
there by the continuity of the solution. We use coordinate freedom to set S|„=o = 0. 

Assumption that the collapse is continuously self-similar, i.e. that there exists a vec- 
tor field C, such that £^gfj,u = '^^g^u, allows the Einstein-scalar equations (^ to be solved 
analytically. The critical solution is given by 



So 


= 0, 




Rq 


= \Ju^ 




$0 


=il„ 

2 


v' 

1 - - 

u. 



(4) 



and is a member of a one-parameter family of self-similar solutions 



s = 


0, 


R = 




$ = 


lln 




2 



uv — pv^, 

2pv + u{l + ^1 + Ap) 



(5) 
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Above fip) is constant component of a scalar field, chosen so that $|„= 
solution converges to the critical solution as p goes to zero: 



M = "2 In 



--Inp + p + 0{p^). 



0, and the 



(6) 



For values p > of control parameter, a black hole is formed; for p < the field disperses, 
and the spacetime is flat again in region u > 0. (Our choice of sign of p is opposite to the 
one in [T^. Also note that the constant field term is omitted there.) 

The existence and position of the apparent horizon is given by condition A{u, v) = 0, 
where 



A{u,v) 



u 



Ry. 



(7) 



Due to the spherical symmetry of the problem it is possible to introduce a local mass M(u, v) 
by 



1 - 



2M 



2g'^'"R^uR,v 



(8) 



The mass contained inside the apparent horizon is the mass of the black hole 



^BH = --Rah, 



(9) 



where i?AH is the radius of apparent horizon given by equation (|^). 



III. LINEAR PERTURBATIONS 



We now apply the perturbation formalism . It is convenient to analyze perturbations 
to the self-similar solution in new variables natural to the problem. For this purpose we 
introduce a "spatial" coordinate x and a "scaling variable" s 



X 



with the inverse transformation 



In 



s = — ln(— u). 



u = —e 



-SI 'Ix 



e^" - 11 



(10) 



(11) 



The signs are chosen to make the arguments of the logarithm positive in the region of interest 
(f > 0, M < 0), where the field evolution occurs. 
In these coordinates the metric (0) becomes 



1 - e-''')ds' - 2dsdx + R' dn\ 



(12) 



and the critical solution (HI) is simply 



0, Ro = 



$0 = X. 



(13) 
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We now consider small perturbations of critical solution (^3]). Quite generally, the per- 
turbation modes can be written as 



S = So+J9^T(x)e'=^ 
i^ = i^o + pr(x)e^-'e^^ 

<^ = %+p<P{x)e^\ (14) 

which amounts to doing Laplace transformation on general perturbations dependent on both 
X and s. It is understood that there can be many perturbation modes, with distinct a, r 
and for each (possibly complex) k in the perturbation spectrum. To recover a general 
perturbation we must take the sum of these modes. Modes with Re > (called relevant 
modes) grow and eventually lead to black hole formation, while modes with Re < decay 
and are irrelevant. 

Let us discuss how the mass-scaling exponent 13 is related to the perturbation spectrum. 
The function A, defined by equation (^, takes the form 

A{x, s, p) = Ao{x) - pa{x)e^\ (15) 

where Aq is the self-similar background term. Recall that A vanishes at the apparent hori- 
zon. Even thought there is no apparent horizon in critical solution, i.e. Aq{x) > 0, the 
exponentially growing perturbation will eventually make A zero, leading to black hole for- 
mation. The mode with the largest eigenvalue, k = max{Re k}, will clearly dominate, giving 
the position of the apparent horizon -Rah(3;ah) by 

-4o(xah) = pa(xAH) 7^ — --• (16) 

Here we made use of the fact that R = e^~^ to the zeroth order in p. Therefore the 
dependence of the black hole mass on p, taken to be the control parameter, is 

Mbh oc i?AH cx (17) 

This is precisely the mass-scaling relation (|l|) with exponent (3 = I/k. For a more detailed 
discussion of this approach and its validity see Ref. 

Perturbing the Einstein-scalar equations (^ by a mode (|T^ , we obtain a system of linear 
partial differential equations. After a change of variables {u,v) to {x,s), the s-dependence 
separates, and we are left with a following system of linear ordinary differential equations 
for cr, r, and in the independent variable x: 

(1 - e-2^)r" + 2{k - e-^^y + A{k - l)r + 4a = 0, (18a) 
(1 - e-2">" + 2{k + e-2">' - 4a + 2(1 - e"'")./.' + 2k<j) + 2e-2V + 2(2 - k)r = 0, (18b) 

r" - 2a' + 20' = 0, (18c) 
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+ 2(e2^ - 2 + e~2x)0/ _ 4^(1 _ ^2x^0 ^ (^;Lg^) 



(1 - e-^^)0" + 2A;0' + 2k(f) + 2(1 - e"^^)/ + 2A;r = 0. (18e) 

Here a prime denotes the derivative with respect to x. 

Boundary conditions for the system (|T8]) are specified at a; = and x = 00. We require 
that perturbations grow slower than background at large x (i.e. that a and r are bounded, 
and (f) grows slower than x) for the perturbation expansion to be valid. Because the solution 
must be continuous at the hypersurface f = 0, we have r(0) = 0(0) = 0. Since we used 
coordinate freedom in u to set S to zero on that surface, we also have cr(0) = 0. Thus, the 
boundary conditions are 

a{0) = r(0) = m = 0, 
cr(x), r(x), 0(x) grow slowly as x — >• 00. (19) 

Equations (|18D together with boundary conditions (|19D constitute our eigenvalue problem. 

The infinitesimal coordinate transformation v ^-^ v — 2({v) corresponding to the remain- 
ing coordinate freedom in v will give rise to unphysical gauge modes 

*E = C(v) = Cle-'(e'" - 1)1, 

5$ = = -e^-'"C[e"'(e'" - 1)1. (20) 

(u-v) -^i y n \ J 

Since on a f = hypersurface 5S, 5R, and 5$ must vanish identically, C(0) = C'(0) = 0, so 
a Taylor expansion of C(f) around zero starts from v"^ term. Therefore these gauge modes 
correspond to negative eigenvalues fc, and are irrelevant. 



IV. SOLUTION OF PERTURBATION EQUATIONS 

Observe that there is an obvious solution of perturbation equations (|18|), which can be 
obtained from self-similar solution (|) by expanding around p = 

S = 0, 



R = \/u^ — uv (1 



^ I = e"-' fl - 2psinh2 X 

2 — uv J ^ 



$ = - In 
2 



1-^ 



p V — 2uv 
2 — uv 



x + psinh2x. (21) 



Clearly, the self-similar mode has eigenvalue = 0, hence any relevant mode will dominate 
the calculation of mass-scaling exponent. 
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The perturbation equations (|T8D inherit their structure from general Einstein-scalar equa- 
tions (H); equations (|18ai|18b|J18e| ) are dynamical equations, while equations (|18c| , |18d| ) are 
constraints. Equation ( |18b| ) is redundant 



(186) = [(^y - 2(|18^ - (1 - e-2")(|Tp' - 2(e-2" + k- 1)(|T8^) + 2(|8iD]/2, (22) 



so it will be discarded. The simple integrable form of equation ( |18c| ) allows elimination of 
one of the unknowns, say a, from the other equations by 



a = r'/2 + + C. 
Thus we are left with two second order ODE's for r and 0: 

(1 - e-2^)r" + 2{k + l- e-2^>' + 4{k - l)r + 40 + 4C = 0, 



(23) 



(24a) 



:i-e 



-2x\ /// 



+ 2k(f)' + 2k(f) + 2(1 - e^^'^y + 2kr = 



(24b) 



with boundary conditions (0), and one constraint (|18d| ), which we equivalently rewrite as 



first order ODE with the use of equation ( 18a ) 



{k - 2)(1 - e-2^)r' + 2{k'^ - fce'^^ - 2k + 2)r 
+ 2(1 - e-2^)0' + A{k - 1)0 + 2C{k + fce'^^ - 2) = 0. 



(24c) 



Equations (|24a|J24b|) can be separated with introduction of the auxiliary function z = r — • 
reducing to 



1 - e-2^)^" + 2kz' + 2{k - 2)z = -AC, 



(25a) 



(1 - e-2^)r" + 2(A; + 1 - e'^^)/ + Akr = A{z - C). 
We convert this into algebraic form using change of variable 



2x ' 1 

y = e , x = -lny. 



Finally, equations (p5D and constraint ( p4cD become 



1/(1 - y) + [1 - (A: + l)y] ^ - {k/2 - l)z = C, 



dz 
dy 



(25b) 



(26) 



(27a) 



d^T dv 
1/(1 - y) ^ + [2 - (A; + 2)y] — - kr = C - z, 



(27b) 



dz dv 
2?/(l -y)-^- Mk - l)z - ky{l -y) — -k{l- ky)r + C[k + {k - 2)y] = 0, (27c) 



and boundary conditions (|T9D are specified at y = 1 and y = oo by 

z(l)=r(l) = 0, ^(1) + C = 0, 
dy 

z{y), r{y) grow slowly as y ^ oo. (28) 

"Grow slowly" means that r{y) must be bounded, and z{y) grows at most logarithmically, 
with ydr/dy + r — z bounded, as y — >■ cxd. The last condition follows from a being bounded 
at infinity. 

Equation (|27c|) is indeed a constraint, since rff pTtj) / dy = (|27q]) — fc (|27fc| ), and is automat- 
ically satisfied for all y if 

{k-l)[kr{l)-2z{l) + 2C] = (29) 
is satisfied initially at y = 1. Imposing initial conditions (^), the constraint ( p9D yields 

C{k - 1) = 0, (30) 

i.e. C = unless k = 1. 

Observe that equations (P7a|) and ( P7b|) are inhomogeneous hypergeometric equations. 
Their properties are well-known and described in a number of books on differential equations 
(see, for example, |13|), so the system (pTj) can be analytically solved and the perturbation 
spectrum exactly determined. 

Equation ( p7a| ) is the homogeneous (except when k = 1) hypergeometric equation 

d'^ z dz 

with coefficients 

c = 1, a + b = k, ab = k/2 — l, 



a= 1/2(A;- VA;2-2fc + 4), b= 1/2 {k + Vk^ - 2k + A). (32) 

It has singular points at y = 0, 1, oo, and its general solution is a linear combination of any 
two different solutions from the set 

zi = F{a,b; k;l- y), 

Z2 = {1- yf~^F{l - a, 1 - 6; 2 - A;; 1 - y), 

= {-yy^F{a, a;a + l-b; y~^), 
z, = i-y)-''F{b,b;b+l-a;y-'), (33) 

where F{a, b; c; y) is the hypergeometric function; F[a, b; c; 0) = 1. Any three of the functions 
(0) are linearly dependent with constant coefficients, for example 

- r(2-fc)r(6-a) r(2-A:)r(a-5) 
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Functions zi, Z2 and z^, Z4 have the asymptotic behavior 

zi ~ 1, Z2 — {I — vY^^ near y = 1, 
Z3 ^ (-y)"", ^ {-yy'' near ?/ = 00. (35) 

The case when = 1 is special. In this case the solution (^) degenerates and zi = Z2 
identically. In order to deal with this situation we continue to use Zi as earlier, but denote 
by Z2 a solution independent of Zi. It is easy to verify that Z2 is logarithmically divergent 



at y = 1. Note also that for /c = 1 we can have C 7^ 0, so that the equation (|27a]) is 
inhomogeneous, and we must add the particular solution z = C to the general solution 
above. 

Observe that at y = 1, Z2 diverges as a power for Refc > 1, logarithmically for k = 1, 
does not have a limit for ReA; = 1, ImA; 7^ 0, and converges (to zero) only for Refc < 1. 
The solutions that satisfy the initial condition z{l) = are z = CZ2, Kek < 1 and z = 



C{1 — zi), k = 1. From equation (|3^) we see that Z2 (and zi) is connected to zs, Z4 by 
a linear relation with non-zero coefficients, so the boundary conditions at infinity will only 
be satisfied if both zs, z^ do not blow up, i.e. if Re a > 0, Re 6 > 0. The curve Rea = 0, 
written in terms the real and imaginary parts of k, has the form 

2 Re fc (2 - Re k) 
^'"^^^ = l-(2Refc)-i ' 

which divides the complex fc-plane into two regions shown in Fig. 0; to the left of the curve 
Rea < 0, and to the right Rea > 0. So to satisfy all boundary conditions (28), k must lie 
in the shaded region K of the complex plane, and for any k & K, 

z = {l- yY^'Fil -a,l-b;2-k;l-y) (37) 

is a solution of equation ( p7a|) with boundary conditions (|28|) . 
Equation (ETO) is also hypergeometric, with coefficients 



c = 2, a = 1, b = k. (38) 
The general solution of the homogeneous equation is 

n=y-\ r2 = y-\l-yY-\ (39) 
and the solution of inhomogeneous equation is easily constructed from (p9D by 



y{l-k) 



y 



1 



r,\i-k 



(40) 



The limits of integration are chosen so that the boundary conditions at ?/ = 1 are satisfied. 
Asymptotic (^) and equation (^OD give the behavior 

y~\l-yf-' ^ 
r ~ near y = I. (41) 

Clearly, r(l) = dr/dy{l) = 0. Also, r is bounded at infinity since z is, for the constraint 
( P7c| ) requires kdr/dy{oo) — 2dz/dy{oo) = 0. Therefore, (p^) and ( ^OD for A; G -ft' is a 
solution of perturbation equations (pTl) with boundary conditions (^), and the region K is 
the perturbation spectrum, which turns out to be continuous. 
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V. CONCLUSION 



We have perturbed the continuously self-similar critical solution of the gravitational 
collapse of a massless scalar field (Roberts solution), and solved the perturbation equations 
exactly. The perturbation spectrum was found to be not discrete, but continuous, and 
occupying the region of the complex plane 



- < Re A; < 1, \lmk\ > 



Re (2 - Re k) 
^ 1- (2Refc)-i 



(42) 



which, according to 0, suggests non- universality of critical behavior for different ingoing 
wavepackets. The complex oscillatory modes might lead to decay of a continuously self- 
similar solution (^) into discrete self-similar choptuon observed in 

The eigenvalues k approach sup Re A; = 1, which corresponds to the mass-scaling ex- 
ponent /5 = 1. This is different from the exponent /3 = 1/2 found in |[T0|-[12| using the 
self-similar solution. 
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FIGURES 



Im k 
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FIG. 1. Complex perturbation spectrum. Values of k to the left of the solid line are prohibited 
by the boundary conditions at infinity, to the right of the broken line by the initial conditions at 
y = 1. Values in the region of intersection (the shaded region K) are allowed, and constitute the 
perturbation spectrum. 
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